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Ultra-Luminous X-ray sources (ULX) are X-ray binaries with L x > 10 39 ergs 
s _1 . The most spectacular examples of ULX occur in starburst galaxies and are 
now understood to be young, luminous High Mass X-ray Binaries. The condi- 
tions under which ULX form are poorly understood, but recent evidence suggests 
they may be more common in low metallicity systems. Here we investigate the 
hypothesis that ULX form preferentially in low metallicity galaxies by search- 
ing for ULX in a sample of Extremely Metal Poor Galaxies (XMPG) observed 
with the Chandra X-ray Observatory (CXO). XMPG are defined as galaxies with 
log(0/H) + 12 < 7.65, or less than 5% solar. These are the most metal-deficient 
galaxies known, and a logical place to find ULX if they favor metal poor systems. 
We compare the number of ULX (corrected for background contamination) per 
unit of star formation (N ULX (SFR)) in the XMPG sample with ~N ULX (SFR) in 
a comparison sample of galaxies with higher metallicities taken from the Spitzer 
Infrared Galaxy Sample (SINGS). We find that ULX occur preferentially in the 
metal poor sample with a formal statistical significance of 2.3a. We do not see 
strong evidence for a trend in the formation of ULX in the high metallicity sam- 
ple: above 12+log(O/H)~8.0 the efficiency of ULX production appears to be flat. 
The effect we see is strongest in the lowest metallicity bin. We discuss briefly the 
implications of these results for the formation of black holes in low metallicity 
gas. 



1. Introduction 



Ultra-Luminous X-ray sources (ULX) are X-ray binaries with L x > 10 39 ergs s _1 most 
commonly found in star forming and starburst galaxies. These sources have attracted con- 
siderable attention in recent years because they have broad-band X-ray luminosities many 



times the Eddington limit for a neutron star or stellar mass black hole (Miller & Colbert 



2004 Roberts |2007 |Feng &: Soria|20l"T ). Some of these extreme objects may be intermediate 



mass black holes (M>5OOM , (Colbert k Mushotzky 1999 Farrell et al. 2009 Sutton et al. 



2012) . It seems likely, however that most ULX are stellar X-ray sources which are either 



radiating in excess of the Eddington limit and/or have black-hole masses somewhat higher 



than is commonly seen in black hole candidates in Milky Way binaries (M>1O , Roberts 



(2007); Gladstone at al. (2009); Zampieri & Roberts (2009) ). ULX are found in a wide va- 
riety of systems - spirals, interacting starbursts and dwarf galaxies. However, the conditions 
under which ULX form are poorly understood. 



It is well established that the number of ULX in star forming galaxies scales with the 
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star formation rate (Grimm et al. 2003 Ranalli, Comastri, & Setti 2003 Mapelli et al. 



2010| |Mineo et al. |[20l2| . In addition, there are are several lines of evidence to suggest that 



ULX form preferentially in low metallicity gas. Kaaret, Schmitt, & Gorski (2011) found that 



the ratio of X-ray luminosity to star formation rate is an order of magnitude larger in low 
metallicity Blue Compact Dwarf Galaxies than for solar metallicity star forming galaxies. 



Swartz et al. (2008) found that the occurrence rate per unit galaxy mass is higher in dwarfs 



than in more massive galaxies (see also Walton et al. (2011)). This surprising correlation is 
explained if ULX favor the metal poor environments found in dwarfs. Evidence that there 



is a direct connection between metallicity and ULX production was presented by Mapelli et 



al. (2010) who find an anti-correlation between the number of ULX and metallicity based on 



a sample of 64 galaxies. A further study by Mapelli et al. (2011) with the addition of two 



XMPG came to the same conclusion. Spectroscopy of ULX counterparts and surrounding 



gas also suggest that ULX are formed from stars in metal poor gas (e.g. Soria et al. (2005); 



Liu et al. (2007)). 



Here we investigate the hypothesis that ULX form preferentially in low metallicity 
galaxies by searching for ULX in a sample of Extremely Metal Poor Galaxies (XMPG) 
observed with the Chandra X-ray Observatory (CXO). XMPG are defined as galaxies with 
\og(0/H) + 12 < 7.65, or less than 5% solar (Papaderos et al. 2008). These are the most 



metal-deficient galaxies known, and a logical place to find ULX if they favor metal poor 
systems. Our goal is to compare the number of ULX (normalized to the star formation rate 
and accounting for cosmic background sources) in a sample of XMPG with a comparison 
sample of galaxies with higher metallicities taken from the Spitzer Infrared Galaxy Sample 
(SINGS). For a sample of galaxies we define Njjlx(SFR) to be: 



N ULX (SFR) 



T,Nulx-ENbkg 
Y^SFR 



Here YI-^ulx is the total number of ULX found in a sample of galaxies (e.g. the 
XMPG sample), ^2 Nbkg the total number of expected cosmic background sources (with an 
apparent luminosity > 10 39 ergs s _1 ) in the same sample and SFR is the integrated star 
formation rate (in M Q per year) of the sample galaxies. 



2. Galaxy sample 

Our sample is listed in Table [I] It consists of 25 nearby (d<50 Mpc) extremely metal 
poor galaxies ((O/H) + 12 < 7.65). These are the most metal-deficient galaxies known. 
Most of those are Blue Compact Dwarf galaxies. Three galaxies in the sample (I Zwl8, SBS 
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0335-052, SBS 0335-052W) were observed with Chandra in 2000 QThuan et al~l[2004"l ) The 
remainder were observed as part of a Cycle 11 Large Project. The exposure times were set 
to obtain a 3.5<r detection of a point source of luminosity 7.8 xlO 38 ergs cm -2 s _1 



survey is complete down to 7.8 xlO 38 ergs cm 2 s 1 . 
Section |3~T Observational details are given in Table [2j 



i.e. our 

Completeness is further discussed in 



3. X-Ray Obervations 



All observations were obtained with the back-illuminated chip ACIS-S3 except for the 
galaxies SBS 0335-052 and SBS 0335-052W, which were observed with the front-illuminated 
ACIS-I3 camera. The observations were performed in VFAINT mode. We reprocessed raw 
data (level 1 event files) using Chandra Interactive Analysis of Observations software (CIAO), 
version 4.2 and the CALibration DataBase (CALDB), version 4.3.0. Standard routines were 
used to correct for bad pixels, charge transfer inefficiency and time dependent gain. A new 
Level 2 events file was created by filtering for standard grades (0, 2, 3, 4, 6) and rejecting 
grades associated with bad pixels. 

The positions of X-ray sources were determined with the CIAO tool Wavdetect. This 



tool is a wavelet-based source detection algorithm (Freeman et al. 2002). The X-ray image 



is convolved with a wavelet function to produce a "correlation image". Clumps of counts 
(sources) are identified as a local maximum in the correlation image if the scale of the wavelet 
is approximately equal to (or greater than) the size of the clump. Hence, Wavdetect also 
gives an estimate of the size of the source. It is typically run with wavelets of differing scales 
to better detect extended emission. 

Wavdetect was run on the level 2 files in the 0.3-8keV band using scales of 1.0, 2.0, 
4.0, 8.0 and 16.0 pixels. We set a threshold significance for identifying a pixel as belonging 
to the source at 10~ 6 . All the detected sources have scales consistent with their being point 
sources. Sources within the D 2 5 region (defined as the elliptical contour best corresponding 
to the 25 mag arcsec -2 blue isophote de Vaucouleurs et al. (1991)) were considered to be 



associated with the galaxy. The results of the source detection are show in Table [3j 

Source regions were defined based on the Wavdetect source extent. We also create the 
corresponding background regions by setting apertures with larger radii than the sources. 
We perform photometry for these sources, using the dmextract tool. The results of the 
photometry are shown at Table [5j The source counts per pixel follow the Poisson distribution 
and the count error, since the number of counts is very low, the Gehrels approximation: error 
= 1 + y/(N + 0.75) (|Gehrels||l~98~6~|). 
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Table 1: The Extremely Metal Poor Galaxy sample. 



Galaxy 


R.A. 


Dec. 


Distance 


D 25 




12+log(0/H) 




(J2000) 


(J2000) 


(Mpc) 


(arcmin) 


(10 20 cm- 2 ) 




UGC 772 


08h25m55.5s 


+35d32m32s 


3.10 


1.20x0.90 


11.5 


7.24+0.03 


SDSS J210455.31-003522.2 


21h04m55.3s 


-00d35m22s 


13.7 




6.47 


7.26+0.03 


SBS 1129+576 


Ilh32m02.5s 


+57d22m46s 


26.3 


0.75x0.10 


0.87 


7.41+0.07 


HS 0822+3542 


08d25m55.5s 


+35d32m32s 


12.7 


0.27x0.12 


4.82 


7.35 


SDSS J120122.32+021108.5 


12h01m22.3s 


+02dllm08s 


18.4 


0.35x0.12 


1.88 


7.55+0.03 


RC2 A1116+51 


Ilhl9m34.3s 


+51d30ml2s 


20.8 


0.24x0.17 


1.19 


7.51+0.04 


SBS 0940+544 


09h44ml6.6s 


+54dllm34s 


24.7 


1.22x1.54 


1.34 


7.48 


KUG 1013+381 


10hl6m24.5s 


+37d54m46s 


19.6 


0.38x0.25 


1.41 


7.58 


SBS 1415+437 


14hl7m01.4s 


+43d30m05s 


10.4 


0.75x0.15 


1.21 


7.60 


6dF J0405204-364859 


04h05m20.3s 


-36d49m01s 


11.0 


0.48x0.35 


0.88 


7.34 


SDSS J141454.13-020822.9 


14hl4m54.1s 


-02d08m23s 


24.6 


0.34x0.24 


4.17 


7.32 


SDSS J223036.79-000636.9 


22h30m36.8s 


-00d06m37s 


18.0 


0.24x0.20 


5.20 


7.64 


UGCA 292 


12h38m40.0s 


+32d46m01s 


3.5 


1.00x0.70 


1.34 


7.27+0.08 


HS 1442+4250 


14h44ml2.8s 


+42d37m44s 


10.5 


1.13x0.26 


1.53 


7.63 


KUG 0201-103 


02h04m25.6s 


-10d09m35s 


22.7 


0.46x0.19 


2.08 


7.56 


SDSS J081239.52+483645.3 


08hl2m39.5s 


+48d36m45s 


9.04 


0.46x0.23 


4.58 


7.16 


SDSS J085946.92+392305.6 


08h59m46.9s 


+39d23m06s 


10.9 


0.39x0.27 


2.44 


7.45 


KUG 0743+513 


07h47m32.1s 


+51dllm28s 


8.6 


0.70x0.30 


5.17 


7.68 


KUG 0937+298 


09h40ml2.8s 


+29d35m30s 


11.2 


0.62x0.23 


1.87 


7.45 


KUG 0942+551 


09h46m22.8s 


+54d52m08s 


24.4 


0.37x0.18 


1.23 


7.66 


SBS 1102+606 


Ilh05m53.7s 


+60d22m29s 


19.9 


0.60x0.28 


0.59 


7.64+0.04 


RC2 A1228+12 


12h30m48.5s 


+12d02m42s 


21.2 


0.29x0.16 


2.47 


7.64 


I Zw 18 


09h34m02.0s 


+55dl4m28s 


17.1 


0.30x0.20 


1.99 


7.17 


SBS 0335-052 


03h37m44.0s 


-05d02m40s 


52.6 


0.23x0.20 


4.98 


7.25+0.05 


SBS 0335-052W 


03h37m38.4xs 


-05d02m37s 


52.2 


0.10x0.13 


4.96 


7.10+0.08 



NOTES: Distances are taken from the NASA/IPAC Extragalactic Database (NED). D 25 shows 
the diameter of major x minor axis. The Galactic Njj is calculated with Colden tool, 
|http: / / cxc.harvard.edu / too lkit / cold en.jsp ] Metallicities are taken from literature. For the galaxy SDSS 
J210455.31-003522.2 we used the optical diameter from NASA/SAO Image Archive. 



Table 2: X-ray observations 



Galaxy 


Obs.ID 


Date 


Exposure (sec) 


Instrument 


UGG 772 


1 1281 

X XZj(J X 


30-08-2000 


5081 


AGTS-S3 

1 \ V 1 k J UU 


SDSS 1210455 31-003522 2 

UXVUU JZilU^U(J.lJl V VU c/ZjZj . Zj 


11282 


04-00-2000 

Vt: ZUu J 


5007 


AGIS-S3 


SBS 1129+576 

kJ XJ U X X Zj i7 | U 1 \J 


11283 


06-07-2010 

V V V 1 Zj V X V 


1 4755 

Xt: 1 UU 


AGIS-S3 


HS 0822-1-3542 

_L 1 k.J uUZZ | U U t: Z. 


11284 


20-12-2000 

Zjw 1Z ZWU J 


51200 

KJ 1ZUU 


AGIS-S3 


SDSS 1120122 32-1-021 108 5 

UXV U kJ J 1ZU1ZZ . U Zj I VZj 1 1UU. U 


11286 


23-1 1-2000 

ZjU XX ZwUl/ 


8007 


AGIS-S3 


RG2 A1116+51 


1 1 287 


D7_1 1-2000 

Wl XX ZUU <-/ 


1 1 640 

X X Ut:W 


AGTS-S3 

i IV 1 kJ UU 


SBS 0040+544 

UXJU Ut/iU 1 Ut: i 


11288 


18-01-2010 

1(J VJ X iUlU 


16828 


AGIS-S3 

ilVJXJ kJ U 


KUG 1013+381 

1 V \J V.I lulO | UVJ X 


11289 


24-01-2010 

Zt: VJ X iUlU 


0402 


AGIS-S3 


SBS 1415+437 

UXJU 1t:1 iJ 1 t:U 1 


11291 


30-102000 


5114 


AGIS-S3 


6dF 10405204-364850 


11292 


28-05-201 


5010 

UU X V 


AGIS-S3 


SDSS 1141454 13-020822 Q 

UXV U kJ J IIIUI ilu UZUUZZ. J 


1 1 293 

X X Zj Ukj 


l8_12-2000 

XO X Zj ZUW iV 


1 6680 

1UUUW 


ACTS-S3 

1 \ V 1 k J UU 


SDSS 1223036 70-000636 


11294 


25-00-2000 

Zj U WiV ^Uu J 


7715 


AGIS-S3 


UGCA 202 

U VJWil Z<JZ 


11295 


06-1 1-2000 

V/VJ XX ZUu J 


5007 


AGIS-S3 


HS 1442+4250 

XXU X X J_Zj 1 TiZj U w 


11296 


26-1 1-2000 


5188 

' J X VJ VJ 


AGIS-S3 

ilVlU UU 


KUG 0201-103 


1 1297 

X 1Z U I 


06-00-2000 


13500 

X U v 7 • J\ I 


AGIS-S3 

ilVlU UU 


SDSS 1081230 52+483645 3 

kJ XV kJ kJ JUUliUJ . U Z | tiUuUtiU . U 


11298 


18-12-2000 

XO 1Z iuu J 


4777 


AGIS-S3 

i 1 V. 1 kJ UU 


SDSS J085946.92+392305.6 


11299 


18-12-2009 


4782 


ACIS-S3 


KUG 0743+513 


11300 


18-12-2009 


5073 


ACIS-S3 


KUG 0937+298 


11301 


16-01-2010 


5007 


ACIS-S3 


KUG 0942+551 


11302 


19-01-2010 


16020 


ACIS-S3 


SBS 1102+606 


11285 


23-08-2010 


10340 


ACIS-S3 


RC2 A1228+12 


11290 


26-07-2010 


12200 


ACIS-S3 


I Zw 18 


805 


08-02-2000 


40750 


ACIS-S3 


SBS 0335-052 


796 


07-09-2000 


59742 


ACIS-I3 


SBS 0335-052W 


796 


07-09-2000 


59742 


ACIS-I3 



The flux of an X-ray source is proportional to the net count rate, where the constant of 
proportionality depends on the response of the detector and the assumed source spectrum: 

A — — - — , U erg /sec/ cm 2 ) /(counts /sec)} (2) 

CountRate 

We extracted standard spectral responses (the Redistribution Matrix (RMF) and Area 



Response Matrix) and simulated a spectrum for each source. Using the results of Swartz et 



al. (2004), we adopt an intrinsic source spectrum with photon index T=1.7 and assume values 
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Table 3: X-ray source detection. 



Galaxy Number of sources Position 

R.A. (J2000) Dec. (J2000) 

SBS 1129+576 1 Ilh32m02.57s +57d22m36.87s 

RC2 A1116+51 1 Ilhl9m34.13s +51d30ml2.48s 

SBS 0940+544 1 09h44ml6.44s +54dllm34.39s 

I ZW 18 1 09h34m01.97s +55dl4m28.23s 

SBS 0335-052 1 03h37m44.04s -05d02m39.89s 

SBS 0335-052W 1 03h37m38.44s -05d02m37.14s 



for galactic absorption taken from Dickey & Lockman ( 1990 ) . We then estimate the constant 



A from the ratio of the number of counts to the calculated flux in the simulated spectrum. 
Fluxes and luminosities were constructed in the 0.3-8 keV band, using the constant A derived 
from the simulated spectra. They are shown in Table |5j One source, IZwl8, had enough 



counts to perform a spectral fit. The best fit is consistent with the analysis of Thuan et al. 



(2004) and is shown in Table |4j The fluxes and luminosities of the other sources are shown 



in Table i 
Table 4: Spectral fit to I ZW 18. 



Galaxy 


Source count-rate T 
(counts/sec) 


(10 20 cm" 2 ) 


X 2 /dof 


I ZW 18 


0.0111 1.88±g:io 


'• 04 -6.17 


12.69/19 



Table 5: X-ray Photometry 



Galaxy 



Net counts Background counts 



N 



H 



Flux 



SBS 1129+576 
RC2 A1116+51 
SBS 0940+544 
I ZW 18 
SBS 0335-052 



21.47+5.71 
89.04+10.48 
39.20+7.32 
494.44+23.25 
28.21+6.38 



SBS 0335-052W 116.29+11.82 



0.19+1.97 
0.21+1.98 
0.20+1.97 
1.56+2.51 
1.30+2.43 
1.30+2.43 



:io 



20 F cm -2) 

0.87 
1.19 
1.34 
7.54 
4.98 
4.96 



(10 15 erg/cm 2 /sec) (10 39 erg/sec) 
1.14+0.3 0.72+0.19 



5.95+0.70 
1.81+0.34 
7.46+0.35 
0.56+0.13 
2.31+0.23 



2.93+0.35 
1.26+0.24 
4.92+0.25 
1.86+0.42 
7.53+0.77 
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3.1. Completeness and Background Sources 

In this section we demonstrate that our sample of ULX in the metal-poor galaxies is 
complete: i.e. that we are not "missing" ULX due to inadequate exposure. In addition, we 
estimate the number of ULX which are chance coincidences: background AGN that happen 
to align with the optical galaxy. 

The completeness limit is expressed by the source detection probability of a galaxy . 
It is a function of the source and background intensity, measured in counts. The detection 
probability as a function of source counts and background counts/pixel can be parameterized 



by the following function (Zezas et al. 2007): 



A(C) = 1.0 - X C- Xl e- X2C (3) 

where C is the source intensity in counts, Ao, Ai, A2 parameters that depend on the back- 
ground counts per pixel. 

All galaxies of our sample were measured to have background level below 0.025 counts 
per pixel. The best-fit parameters for this background give a detection probability of the 
form: 



A(C) = 1.0 - 11.12C- a83 e- a43C . (4) 

Therefore, by solving the equation for 90% and 50% completeness, the corresponding 
source counts are: C90 = 7.2 counts and C50 = 4.4 counts. 



Table [6] shows that all luminosities of 90% completeness and all luminosities 50% com- 
pleteness, are well under the limit of 10 39 ergs s— 1 cm 2 . This means that we have detected 
at least 90% of all existing ULX sources in the galaxy sample. 

The predicted number of background sources for each galaxy was estimated as follows. 
For each galaxy, we calculate the 0.3-8 keV flux for a source of 10 39 ergs s -1 and then 
transform the flux to the 0.5-2 keV band assuming a power law source spectrum with photon 
index T=1.7 and a foreground N# listed in Table [TJ We then used the LogN-LogS curves 



of Giacconi et al. (2001) for the 0.5-2 KeV band to estimate the number of background 
sources per square degree and normalize by the D25 values to get the absolute number of 
background sources predicted to lie within the optical area of the galaxy. These values are 
given in Table [6} We note that SBS 0940+544 has a predicted background value close to 
one. 
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Table 6: Completeness 



galaxy 


Exposure time 


yu /o 


OU /0 


Nbkg 




(sec) 


(lO^ 8 erg/sec) 


(lO'WS er g/ sec ) 




UGC 772 


5081 


1.66 


1.01 


0.014 


SDSS J210455. 31-003522. 2 


5007 


2.39 


1.46 


0.039 


SBS 1129+576 


14755 


2.40 


1.47 


0.034 


HS 0822+3542 


5120 


2.00 


1.23 


0.004 


SDSS J120122. 32+021108. 5 


8097 


2.66 


1.63 


0.010 


RC2 A1116+51 


11640 


2.37 


1.45 


0.012 


SBS 0940+544 


16828 


2.31 


1.41 


0.776 


KUG 1013+381 


9402 


2.60 


1.59 


0.026 


SBS 1415+437 


5114 


1.35 


0.82 


0.010 


6dF J0405204-364859 


5010 


1.54 


0.94 


0.017 


SDSS J141454.13-020822.9 


16680 


2.31 


1.41 


0.034 


SDSS J223036.79-000636.9 


7715 


2.67 


1.63 


0.012 


UGCA 292 


5007 


0.08 


0.05 


0.041 


HS 1442+4250 


5188 


1.35 


0.83 


0.028 


KUG 0201-103 


13590 


2.41 


1.47 


0.031 


SDSS J081239.52+483645.3 


4777 


1.09 


0.67 


0.008 


SDSS J085946.92+392305.6 


4782 


1.58 


0.97 


0.011 


KUG 0743+513 


5073 


0.93 


0.57 


0.014 


KUG 0937+298 


5007 


1.59 


0.97 


0.015 


KUG 0942+551 


16020 


2.36 


1.45 


0.027 


SBS 1102+606 


10340 


2.44 


1.49 


0.048 


RC2 A1228+12 


12200 


2.37 


1.43 


0.014 


I ZW 18 


25956 


0.72 


0.44 


0.013 


SBS 0335-052 


59742 


4.74 


2.89 


0.070 


SBS 0335-052W 


59742 


4.66 


2.85 


0.020 



4. Estimates of the Star Formation Rate in the XMPG Sample 



It is well established that the number of ULX in a galaxy scales with the star formation 
rate QGrimm et aT~||2003| |Ranalli, Comastri, fc Settl||2003| |Mapelli et al^OlO} |Mineo et al~ 



2012). We therefore need reliable estimates of the star formation rate in the XMPG sample 



in order to determine whether ~Njjlx(SFR) is higher in the metal poor sample than in the 
SINGS sample. We use two methods to estimate the star formation rate: the Far Ultraviolet 
(FUV) luminosity from GALEX and the 24/zm luminosity from Spitzer. 



4.1. GALEX Data 



The FUV emission from star forming regions comes directly from young massive O and 



B stars. We use the method from Hunter et al. (2010) to relate the FUV luminosity to the 
star formation rate in dwarf galaxies: 



SFR FXJV (M Q /yr) = 1.27 x 10" 28 L FUV (er^s- 1 ^- 1 ) (5) 

We obtained GALEX FUV images from the GalexView version 1.4.6 catalog. All galax- 
ies in the XMPG sample were detected with the exception of SDSS J223036. 79-000636. 9. 
This galaxy is excluded from the rest of the analysis in this paper. Source counts were 
extracted using the CI AO dmextract routine. This was done with either circular or el- 
liptical source apertures, depending on the morphology of the galaxy in question, as well as 
an annular background aperture. Both of these apertures were centered on the source, with 
the circular apertures having a radius of either 10 or 15 pixels, and the ellipses having a 
semi-major and semi-minor axis of length 28 and 14 pixels or 15 and 7.5 pixels, respectively, 
depending on the size of the source. All background annuli had an area 8 times that of their 
respective source apertures. The fluxes were computed using the relationship /fuv = (1-4 X 
10~ 15 ) x count rat^| Note that these fluxes are in units of erg cmT 2 s _1 A -1 . We converted 
these values to units of erg cm~ 2 s _1 Hz -1 by multiplying them by A 2 /c, where A is the 
effective wavelength of the FUV detector (= 1528 A), and c is the speed of light in units of 
A s -1 . With the fluxes in these new units, we proceeded to calculate the Lfuv values of the 
sources, which were subsequently corrected for extinction. We then calculated the SFRfuv 
estimates of the galaxies using the extinction-corrected luminosities. 



1 http://galexgi. gsfc.nasa.gov/docs/galex/FAQ/counts_background. html 



- 11 - 



Table 7: UV measurements and SFRs 



Galaxy 


E(B-V) (1) 


Count Rate (2) 




SFRfttv^ 




fmag) 


(ct s" 1 ) 


(x 10 26 erg s _1 Hz -1 ) 


flO -3 M n yr" 1 ) 


TTPP 779 


098 
u.uzo 


O.OOjzu.UO 


n riQ-i-o 001 
u.uyznu.uui 


i riQ-i-o 01 
i.uyznu.ui 


ciDcm 79104^^1 00^99 9 


o nfifi 

u.uuu 


i 91 _|_n On 
1 .Z- IJZU. uu 


4Q-I-0 09 
u.^yznu.uz 


« 99-un 90 
u.zzznu.zy 


C1RC1 1 1 9Q4-^7fi 

0J30 lliJTtJiU 


013 


3 48-1-0 90 
o.^o znu.zu 


3 47-I-0 90 


44 09-I-9 ^9 
^t^t. u zjz z. o z 


0899-l-'? c i49 
no uozz-t-oo^z 


047 

U.U^ 1 


i nQ_i_n no 
i .uoznu.uy 


n Q9-i_n n3 
u.oz jzu.uo 


a 1 1 _i_n Q9 

iinzu. oo 


ClDCm T190199 ^9-1-091108 ^ 


094 


i 9fi-l-0 03 
i.zuznu.uo 


n4+0 09 
u.u^nzu.uz 


o H7-I-0 90 
o.u i inu.zu 


RP9 A 1 1 1 fi-l-^l 


01^ 

U.UIO 


a Qfi-i-n 91 


o c o_i_n i o 

Z.OOILU. ±Z 


09 7 n _|_i c,7 


SRS 0040-1-^44 

Q-L>Q Uy^iUn^O^i^i 


01 3 

U.UIO 




± .o i znu. -Lu 


17 41+1 9c; 


KTTfl 1 01 ^-1-^81 


01^ 


7 n Q-|-n 9fi 


^i.Z^iZCU. 1U 


JO.O 1 _T_ l.OI 


qnq i 41 c:_i_/LQ7 


00Q 
u.uuy 


i o «o-i-n 30 
iz.uyznu.ou 


i QQ_|_n 
i .yoinu.uo 


9/L ^^-1-0 


firlP 7040^904 ^48^0 
our j u^uozu^-ou^ooy 


OOfi 
u.uuu 


r 9Q-|-n 99 
o.zymu.zz 


n fifi-l-0 03 
u.oUjzu.uo 


in 1 9-1-0 
iu. izmu.^o 


ciDcm 11414^4 13 090899 Q 


0^8 
u.uoo 


n n Q-un 09 
u.uoznu.uz 


pi 7Q-un n9 
u. i oznu.uz 


j.OOIlU.OU 


TTHPA 9Q9 


01 n 
U.U1U 


i . iohzu. lo 


n 1 3-I-0 003 

U. lOILU.UUO 


1 (^o-un OA 
1 . uZfUzu .u^± 


TTO "1 A A O I /lorn 

Ha 1442+4250 


n m o 

0.013 


n nn i n oo 

9.00+0.22 


1 /i o i n n /i 

1.42+0.04 


"i n 1 n ^ r 

18.08+0.45 


KUG 0201-103 


0.021 


1.30+0.09 


1.02+0.07 


13.00+0.87 


SDSS J081239.52+483645.3 


0.051 


0.86+0.02 


0.13+0.003 


1.71+0.04 


SDSS J085946.92+392305.6 


0.026 


0.71+0.01 


0.14+0.001 


1.72+0.02 


KUG 0743+513 


0.070 


12.58+0.34 


1.66+0.04 


21.12+0.57 


KUG 0937+298 


0.018 


5.21+0.16 


0.97+0.03 


12.34+0.37 


KUG 0942+551 


0.012 


1.48+0.12 


1.26+0.10 


16.01 + 1.27 


SBS 1102+606 


0.006 


4.19+0.17 


2.27+0.09 


28.79+1.14 


RC2 A1228+12 


0.027 


2.69+0.04 


2.03+0.03 


25.69+0.39 


I ZW18 


0.032 


11.52+0.08 


5.60+0.04 


71.12+0.52 


SBS 0335-052 


0.047 


5.52+0.20 


28.41+1.00 


360.86+12.76 


SBS 0335-052W 


0.046 


0.56+0.01 


2.84+0.04 


36.07+0.53 



Notes: ^ Reddening magnitude of the source, used to calculate source extinction; ( 2 ) 
Observed background-subtracted source count rate; ® Extinction-corrected far-UV 
luminosity of the source; ( 4 ) Star-formation rate, calculated using the formula from Hunter 

et al. (2010). 
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4.2. Spitzer Data 



The 24/im emission in galaxies comes from single photon transient heating of small 



grains and can be used as a tracer of recent star formation. Calzetti et al. (2007) derived the 



following relation between the star formation rate (SFR) and 24/im emission from calibrating 
Hjj regions in nearby galaxies: 



SFR 1R (M e /yr) = 1.31 x 10- 38 [L 24/ , m (er^ S / S ec)] a885 , (6) 

where for intermediate luminosity galaxies: 

lxlO 40 < L 24Mm <3xl0 44 erg/sec 

We acquired data from the Spitzer space telescope that were obtained with the MIPS 
instrument in the 24 microns band. We obtained the post-BCD (post Basic Calibrated Data) 
data for the 9 galaxies in the sample that have available MIPS 24/im data. The data sets 
consist of several exposures that are interpolated into one mosaic image. From these images 
we measure the surface brightness in units of MJy/sr. The post-BCD product mosaic pixel 
sizes are 2.42x2.45 arcsec 2 for the 24 microns detector. 

We perform aperture photometry using funtools. We define an aperture for the source 
with an ellipse that includes the total flux of the source and an elliptical annulus for the 
background region as depicted in Figures 15 - 16 for the 9 galaxies. In order to calculate the 
total flux density of each source, we sum the total flux of the source, subtract the background, 
on a set of pixels and multiply by the number of steradian per pixel. We convert the surface 
brightness to the flux density with the following formula: 

f v =Y 2.45 x 2.45^^ x 0.023504 FluX J mJyl (7) 
z — f pixel arcsec 2 

pixels 

The flux density is converted to monochromatic flux: 

F=jf v (X). (8) 



The uncertainty on the flux is calculated from the mosaic variance image by adding in 
quadrature the uncertainty of all pixels within the source aperture. The results are listed at 
Table E 
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4.3. Comparison of Spitzer and GALEX star formation rates 

Several galaxies in the XMPG sample have both GALEX UV images and Spitzer 24 /im 
data. It is therefore instructive to compare SFR derived from the two methods. Figure [T] 
shows the star formation rate derived from the FUV flux with the star formation rate derived 
from the 24/iin flux. There is a clear correlation between the two methods. The rate derived 
from the FUV flux is systematically higher than for the 24/im flux. This is almost certainly 
because 24/zm emission underestimates the star formation rate because of the very low dust 
content of the XMPG. We adopt the GALEX derived star formation rate because (1) the 
Spitzer 24 /im data likely underestimates the star formation rate, and (2) it is available for 
all the galaxies. We note that SBS 0335-052 appears to have very high SFR in comparison 
with the other XMPGs. 



5. The SINGS sample. 



In order to determine whether there is a statistically significant excess of ULX in the 
most metal poor systems, we require a comparison sample of galaxies with well determined 



star formation rates and metallicities. The Spitzer Infrared Nearby Galaxies Survey (Ken- 



nicutt et al. 2003 ) is ideal for this purpose. It is a comprehensive imaging and spectroscopic 



study of 75 nearby galaxies with distances < 30 Mpc. The morphological types of the sample 
range from elliptical to irregular. The SINGS sample does not include absolute extremes in 
properties that can be found in larger volumes such as ULIRGs, luminous AGNs or extremely 
metal poor galaxies. A distance-limited sample of the SINGS survey has been observed by 
Chandra (Jenkins et al in preparation). We use a sub-sample of SINGS galaxies selected 



by Calzetti et al. (2007) which have Hubble Space Telescope (HST) NICMOS images in the 
Paa hydrogen emission line (1.8756 jum), Ha observations and Chandra observations. The 
Ha and the Paa lines are used to measure the extinction correction. Furthermore, the Paa 
line is used to calibrate the mid-infrared emission. There are 33 galaxies in the Calzetti 



et al. (2007) sample, 26 of which also form part of the Chandra survey. The galaxies are 



divided into three groups according to their oxygen abundance: High-metallicity galaxies 
[12+log(0/H) > 8.35], Intermediate-metallicity galaxies [8.00 < 12+log(0/H) < 8.35], and 
Low-metallicity galaxies [12+log(0/H) < 8.00]. For reference, the XMPGs have metallicities 
12+log(0/H) < 7.65. We list the main characteristics of the SINGS sample in Table 19} 



The metallicities given in Table [9] are given as an upper and lower bound and taken from 



Moustakas et al. (2010). The lower bound is the metallicity derived using the calibration 



of Pilyugin & Thuan (2005) (hereafter PT05) and the upper bound using the calibration of 



Kobulnicky & Kewley (2004) (hereafter KK04). As discussed in Moustakas et al. (2010), the 
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Fig. 1. — Comparision of the star formation rate determined from the FUV luminosity and 
24- /im luminosity. The line denotes SFR/R=SFRp;yv'- 
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PT05 calibration is based on empirical abundance measurements of individual HII regions 
and the KK04 calibration on photoionization model calculations. The KK04 calibration 
gives abundances that are systematically higher than PT05. In this paper we choose the 
PT05 calibration since this method was used to derive the XMPG abundances. The results 
of this paper are unchanged if we use KK04 or an average of the two. 

We estimate the star formation rate of the SINGS galaxies using the calibration of 



Calzetti et al. (2010): 



SFR(M Q /yr) = 5.5 ■ 10~ 42 x [L(Ha) obs + 0.02L(2Afim)} 



(9) 



Note that the L(Ha) ft s is the observed Ha luminosity without correction for internal 
dust attenuation. 



We use the measurements of the Ha and 24/zm fluxes of Dale et al. (2007) ) and 



Kennicutt et al. (2008) & Kennicutt et al. (2009). The numbers of ULX for the SINGS 



sample have been provided by Chandra SINGS team (Jenkins et al in preparation). 



6. Comparison between the SINGS and XMPGs samples 

In this section, we investigate the relationship between metallicity and Nulx{SFR). 



Table [10] compares N ULX (SFR) for the XMPG sample and the SINGS sample. The SINGS 
sample is divided into 3 sub-groups according to metallicity values (High, Intermediate and 
Low metallicity). In addition, we show Nulx(SFR) obtained by combining the SINGS low 
metallicity sample and the XMPG sample. The difference in Njjlx{SFR) between the high 
metallicity SINGS galaxies and the low metallicity galaxies (comprising the XMPG and low 
metallicity SINGS galaxies) is 2.3a. The low metallicity sample has a small number of in- 
dividual galaxies with very high N^x/SFR values (these have low SFR and one ULX). A 
Kolmogorov-Smirnov (K-S) test gives the probability that the two distributions (low metal- 
licity and SINGS) come from the same parent population as 0.18. Finally, Figure [2] shows 
Nulx(SFR) as a function of metallicity. The high metallicity SINGS galaxies are plotted 
as individual points, and the XMPG and SINGS low metallictiy galaxies are combined. We 
note that for galaxies with no ULX, the background-subtracted number of ULX is negative, 
which is unphysical. There is a marked increase in Njjlx{SFR) in the low metallicity galax- 
ies. We conclude that ULX form preferentially in low metal systems, with the caveat that 
the formal significance of this result is low. 



We do not find any any evidence for a trend in N ULX (SFR) with metallcity 12+log(O/H)>8.0. 
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Fitting the data points in Figure [2] above 12+log(O/H)=8.0 with both a flat line (slope 0) 
and a straight line with non-zero slope we find that the more complex model (non-zero slope) 
is slightly preferred on the basis of a x 2 fit. We use an F-test to determine whether the model 
with a slope is significantly better than the flat line (the null hypothesis is that a slope does 
not give a statistically better fit). The F-test gives a significance of 0.35, confirming that a 
slope is not required. 



Mapelli et al. (2010) use a larger sample of galaxies and find an anti-correlation between 
Nulx(SFR) and metallicity (their Figure 5 is directly comparable to Figure [2] of this pa- 
per). As discussed in the previous paragraph, there may be a similar trend in the SINGS 
galaxies but the effect is small below 12+log(O/H)=8.0 and the scatter is large. We cannot 
unambiguously confirm this result. The preference for ULX to form at the low metallicities 
is most apparent in the lowest metallicity bin. 



6.1. Summary and Discussion 



In this paper, we present the results of a Chandra survey to search for ULX in the most 
metal poor galaxies. We find that, compared to a comparison sample of high metallicity 
SINGS galaxies, the low metallicity galaxies are more likely to host a ULX. The number of 
ULX normalized to the star formation rate is ~ 0.17 for the high metallicity galaxies and 7.0 
for the low metallicity systems (12+log(O/H)<8.0). The formal significance of this result is 



low. This study broadly agrees with the results of Mapelli et al. (2011) who also find that 



the number of ULX increases in lower metallicity galaxies. Mapelli et al. (2011 ) also claim to 



see a trend in the numbers of ULX as a function of metallicity. As demonstrated in Figure [2j 
we do not see strong evidence that ~Njjlx(SFR) decreases as a function of metallicity above 
~ 8.0. This result suggests there maybe a "threshold" for more efficient ULX formation at 
about 12+log(0/H)<8. 

In a recent paper, we ( Prestwich et al.||2012 ) compared the ULX population of two colli- 
sional ring galaxies, the Cartwheel and NGC 922. The Cartwheel has a relatively low metal- 
licity (12+log(0/H)~8.1) and NGC 922 has a metallicity close to solar (12+log(0/H)~8.81. 
We found that the number of ULXs in NGC 922 and the Cartwheel scales with the star for- 
mation rate: we do not find any evidence for an excess of sources in the Cartwheel relative to 
NGC 922. The N ULX (SFR) values for NGC 922 and the Cartwheel are the same (within the 
errors) as the SINGS galaxies. This is demonstrated in Figure [3] where we plot Nulx(SFR) 
vs. metallicity for the SINGS galaxies, the Cartwheel and NGC 922. The KK04 calibration 
is used for this plot as it is more appropriate for the Cartwheel and NGC 922 measurements. 
Figure [3] adds to the evidence that the high number of ULX in the Cartwheel is due to it's 
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Fig. 2. — Njjlx(SFR) for individual SINGS galaxies, intermediate metallicity galaxies and 
the combined metal poor and XMPG. This plot uses the PT05 metallicity calibration. 
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Fig. 3.— N ULX (SFR) for individual SINGS galaxies, NGC 922 and the Cartwheel. This 
plot uses the KK04 metallicity calibration. 
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very high star formation rate, and not a metallicity effect. 



There are two hypotheses which might explain the excess of ULX in metal-poor galax- 



ies. 



Linden et al. (2010) studied the effect of metallicity on HMXB production and found 



that the number of ULX in low metal systems increased dramatically below a threshold of 
Z/Zq <10%. The high X-ray luminosities derive from Roche Lobe Overflow onto a black 



hole, which typically has a mass M< 1OM . The threshold predicted by Linden et al. (2010) 



is consistent with the increase in Nulx(SFR) we see in galaxies with 12+log(0/H)<8. An 
alternate hypothesis is that higher mass black holes are able to form in lower metallicity 



gas (Mapelli et al. 2009 Zampieri & Roberts 2009) leading to higher X-ray luminosities in 



HMXB. Our results are consistent with both of these scenarios. 
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Table 8 



: 24 microns measurements 


and SFRs 






Galaxy 


Flux density 


Flux 




SFR 




(mJy) 


(10~ 13 erg/sec/cm 2 ) 


(10 40 erg/sec) 


{M yr- V ) 


HS 0822+3542 


2.23+0.09 


2.78+0.12 


0.54+0.02 


0.002+0.001 


SBS 0940+544 


1.60+0.10 


2.01+0.13 


1.47+0.09 


0.005+0.003 


KUG 1013+381 


14.54+0.17 


18.2+0.21 


8.37+0.09 


0.022+0.002 


SBS 1415+437 


16.67+0.15 


20.8+0.19 


2.70+0.02 


0.008+0.001 


UGCA 292 


<0.87±0.31 


1.09+0.34 


0.08+0.03 


0.0003+0.0010 


HS 1442+4250 


3.23+0.11 


4.03+0.14 


0.53+0.02 


0.002+0.001 


SBS 1102+606 


1.33+0.10 


1.67+0.11 


0.79+0.06 


0.003+0.002 


I Zw 18 


4.84+0.11 


6.05+0.14 


2.12+0.05 


0.006+0.001 


SBS 0335-052 


67.25+0.20 


84.10+0.25 


278.94+0.84 


0.480+0.013 


SBS 0335-052W 


0.14+0.04 


0.018+0.005 


0.58+0.17 


0.002+0.005 



NOTES: The SFR errors were calculated with error propagation. For the galaxy UGCA 292 we calculate 
the upper limits. 
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Table 9: Characteristics of the SINGS sample. 



Galaxy 


Morphology 


Distance 


12+log(0/H) 


L(Ha) o6s 


SFR 


Nf/LX 


Nbkg 






(Mpc) 


(10 40 erg/sec) 


(M /yr) 






High-metallicity galaxies 


NGC 0925 


SAB(s)d 


9.12 


8.24-8.78 


6.01 


1.634 


1 


0.566 


NGC 2403 


SAB(s)cd 


3.5 


8.31-8.81 


5.12 


1.225 


1 


0.474 


INGC 2841 


OA/ \1 

bA(rjb 


9.8 


8.52-9.19 


3.39 


3.17 


2 


n no o 

0.282 




oAc 


3.5 


8.30-8.98 


0.87 


n Q 9 A 

0.334 


U 


U.028 


1NGC 3184 


bAB(rsjcd 


11 1 A 

11.10 


8.48-9.14 


7.36 


3.31 


2 


n /^o o 

0.628 


M/^n oi no 

INGC 3198 


br>(rs)c 


1 o /? o 

13.68 


8.32-8.87 


6.85 


3.71 


-I 

1 


n a no 

0.493 


NGC 3351 


SB(r)b 


10.1 


8.60-9.22 


3.71 


3.907 





0.093 


NGC 3627 


SAB(s)b 


8.7 


8.49-9.10 


10.9 


11.362 


2 


0.312 


INljrL/ 4boy 


bAr>(rs Jed. 


ll.l 


8.25-8.79 


10.6 


O K A K 


9 

6 


n k*7Q 
(J.o to 


NGC 4569 


SAB(rs)ab 


16.58 


8.56-9.19 


4.84 


6.731 


1 


1.015 


NGC 4625 


SAB(rs)m 


9.17 


8.27-9.04 


0.61 


0.215 





0.013 


NGC 4736 


(R)SA(r)ab 


5.3 


8.31-9.01 


3.45 


2.709 


3 


0.357 


NGC 4826 


(R)SA(rs)ab 


5.6 


8.59-9.2 


4 4.05 


2.732 





0.209 


NGC 5055 


SA(rs)bc 


7.82 


8.42-9.13 


7.71 


6.172 


2 


0.613 


NGC 5194 


SA(s)bc 


8.2 


8.54-9.18 


20.10 


15.847 


3 


0.569 


NGC 6946 


SAB(rs)cd 


5.0 


8.40-9.04 


16.5 


16.437 


1 


0.4258 


NGC 7331 


SA(s)b 


15.1 


8.40-9.05 


13.4 


15.893 


3 


0.850 






Intermedial 


;e-Metallicity Galaxies 








NGC 1705 


SA0- 


5.1 


8.20-8.43 


0.9 


0.075 





0.009 


IC 2574 


SAB(s)m 


2.8 


7.94-8.26 


0.97 


0.128 





0.084 


NGC 4236 


SB(s)dm 


4.45 


8.07-8.56 


1.3 


0.183 





0.411 


IC 4710 


SB(s)m 


7.8 


8.11-8.62 


1.6 


0.231 


1 


0.071 


NGC 6822 


IB(s)m 


0.47 


8.04-8.67 


0.07 


0.015 





0.013 


Low-Met allicity Galaxies 


Ho II 


Im 


3.5 


7.68-8.07 


0.64 


0.073 


1 


0.087 


DDO 053 


Im 


3.56 


7.77-8.13 


0.08 


0.011 





0.003 


Ho IX 


Im 


3.3 


7.61-7.98 


0.01 


0.009 


1 


0.008 


NGC 5408 


IB(s)m 


4.8 


7.81-8.23 


1.16 


0.023 


1 


0.004 



NOTES: Galaxy morphologies are from the NASA/IPAC Extragalactic Database (NED). Adopted distances as derived by 
|Masters~| ( |2005} . Oxygen abundances are from [Moustakas et al.| | |2010[ ) 
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Table 10: Comparison of the SINGS and XMPGs sample. 





High 


Intermediate 


Low 


XMPGs 


Low+XMPG 




25 


1 


3 


5 


8 




7.5 


0.6 


0.1 


1.3 


1.4 


N ULX (SFR) 


0.17±0.042 


0.65±1.0 


25.0±14.6 


4.5±2.3 


7.0±2.7 



